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ABSTRACT 

We report long-slit spectroscopic observations of the quasar SDSS 
J082303.22-1-052907.6 (2ciV"^3.1875), whose Broad Line Region (BLR) is partly 
eclipsed by a strong damped Lyman-ct (DLA; logA^(HI)=21.7) cloud. This allows us 
to study the Narrow Line Region (NLR) of the quasar and the Lyman-a emission 
from the host galaxy. Using cloudy models that explain the presence of strong 
NV and PV absorption together with the detection of Sill* and Of** absorption 
in the DLA, we show that the density and the distance of the cloud to the quasar 
are in the ranges 180 < nn < 710 cm“^ and 580 > tq >230 pc, respectively. Sizes 
of the neutral(^2-9pc) and highly ionized phases ('^3-80pc) are consistent with the 
partial coverage of the CIV broad line region by the CIV absorption from the DLA 
(covering factor of ~0.85). We show that the residuals are consistent with emission 
from the NLR with CIV/Lyman-a ratios varying from 0 to 0.29 through the profile. 
Remarkably, we detect extended Lyman-a emission up to 25kpc to the North and 
West directions and 15kpc to the South and East. We interpret the emission as the 
superposition of strong emission in the plane of the galaxy up to lOkpc with emission 
in a wind of projected velocity ^500km s“^ which is seen up to 25kpc. The low 
metallicity of the DLA (0.27 solar) argues for at least part of this gas being in-falling 
towards the AGN and possibly being located where accretion from cold streams ends 
up. 

Key words: galaxies: evolution - intergalactic medium - quasars: absorption lines 
- quasars: individual: SDSS J082303.22-k052907.6 


1 INTRODUCTION 

Luminous high-redshift quasars consist of supermassive 
black holes residing at the center of massive galaxies and 
growing through mass accretion of gas in an accretion disk. 
Bright quasars play an important role in shaping their 
host galaxies through the emission of ionizing flux but also 
through launching powerful and high-velocity outflows of 


* Based on data obtained with MagE at the Clay telescope of 
the Las Campanas Observatory (CNTAC Prgm. ID CN2012B-51 
and CN2013A-121) 


gas. These outflows inject energy and material to the disk 
of the galaxy and may be influencing the physical state up 
to larger distances. It has remained unclear however what 
are the mechanisms that drive energy from the very center 
of the active galactic nuclei (AGN) to the outskirts of the 
galaxy. 

Observational evidence for outflows and winds in AGNs 
is seen as prominent radio-jets in radio-loud sources, broad 
absorption lines observed in broad absorption line (BAL) 
quasars, or through the photoionized warm absorber which 
is frequently observed in the soft X-rays (e.g. Crenshaw et al. 
2003). Gravitational micro-lensing studies have shown that 




2 H. Fathivavsari et al. 


the primary X-ray emission region in AGN is of the order of 
a few tens of gravitational radii in size (e.g. Dai et al. 2010) 
and X-ray spectroscopy shows that highly ionized outflows 
launched from this region are seen in high-z quasars (Char- 
tas et al. 2009) and in at least 40% of them with velocities 
up to 0.1 c (Gofford et al. 2013). 

Outflows are observed also on large scales. Mullaney et 
al. (2013) used the SDSS spectroscopic data base to study 
the one-dimensional kinematic properties of [OIII]A5007 by 
performing multicomponent fitting to the optical emission¬ 
line profiles of about 24000, z < 0.4 optically selected AGNs. 
They showed that approximately 17 percent of the AGNs 
have emission-line profiles that indicate their ionized gas 
kinematics are dominated by outflows and a considerably 
larger fraction are likely to host ionized outflows at lower 
levels. Harrison et al. (2014) find high-velocity ionized gas 
(velocity widths of about 600-1500 km s“^) with observed 
spatial extents of (6-16) kpc in a representative sample of 
z < 0.2, luminous (i.e. L[0 lll]>10^^'^ erg s“^) type 2 AGNs. 
Therefore galaxy-wide energetic outflows are not conhned to 
the most extreme star-forming galaxies or radio-luminous 
AGNs. 

If outflows are observed both on small and large scales, 
how the small scale outflows are transported at larger 
distances remains unclear (Faucher-Giguere et al. 2012, 
Ishibashi & Fabian 2015, King & Pounds 2015). This is how¬ 
ever a crucial question as these outflows are massive and 
energetic enough to signihcantly influence star formation in 
the host galaxy and provide signihcant metal enrichment to 
the interstellar and intergalactic media (e.g. Dubois et al. 
2013). At high redshift where quasars are more luminous, 
the consequences of such outflows are of first importance for 
galaxy formation. One way to study the interplay between 
the quasar and its surrounding is to search for Lyman-a 
emission around quasars (e.g. Hu & Cowie 1987, Hu et al. 
1996, Petitjean et al. 1996, Bunker et al. 2003, Christensen 
et al. 2006). These observations reveal gas infalling onto the 
galaxy (Weidinger et al. 2005), positive feed-back from the 
AGN (Rauch et al. 2013) or a correlation between the lumi¬ 
nosity of the extended emission and the ionizing flux from 
the quasar (North et al. 2012). 

Very recently, we searched quasar spectra from the 
SDSS-HI Baryon Oscillation Spectroscopic Survey (BOSS; 
Dawson et al. 2014) for the rare occurrences where a strong 
damped Lyman-a absorber (DLA) blocks the Broad Line 
Region emission (BLR) from the quasar and acts as a nat¬ 
ural coronagraph to reveal narrow Lyman-a emission from 
the host galaxy (Finley et al. 2013; see also Hennawi et al. 
2009). This constitutes a new way to have direct access to 
the quasar host galaxy and possibly, when the size of the 
DLA is small enough, to the very center of the AGN. Out of 
a total of more than 70,000 spectra of z > 2 quasars (Paris 
et al. 2012), we gathered a sample of 57 such quasars and 
followed-up six of them with the slit spectrograph Magellan- 
MagE to search for the very special cases where the DLA 
coronagraph reveals the very center of the host galaxy and 
extended Lyman-a emission. In the course of this follow-up 
program, we found one object SDSS J0823+0529 where the 
DLA does not cover the Lyman-a broad line region entirely 
and reveals the emission of the Lyman-a and C IV narrow 
line regions. We show here that this is a unique opportu¬ 
nity to study the link between the properties of the central 


regions of the AGN to that of the gas in the halo of the 
quasar. 

The paper is organized as follows. In Section 2 we 
describe the observations and data reduction. We derive 
properties of the gas associated with the DLA (metallic- 
ity, ionization state, density, distance to the quasar, typical 
size) in Section 3. We discuss the properties of the quasar 
narrow line region and of the extended Lyman-a emission 
in Sections 4 and 5, respectively. We then finally, present 
our conclusions in Section 6. In this work, we use a stan¬ 
dard CDM cosmology with Da = 0.73, Dm = 0.27, and 
Ho =70 km s“^ Mpc“^ (Komatsu et al. 2011). Therefore 
1 arcsec corresponds to about 7.1 kpc at the redshift of the 
quasar (z = 3.1875 see below). In the following we will use 
solar metallicities from Asplund et al. (2009). 


2 OBSERVATIONS AND DATA REDUCTION 

The spectrum of the quasar SDSS J0823-I-0529 was observed 
with the Magellan Echellete spectrograph (MagE; Marshall 
et al. 2008) mounted on the 6.5 m Glay telescope located at 
Las Campanas Observatory. MagE is a medium-resolution 
long-slit spectrograph that covers the full wavelength range 
of the visible spectrum (3200 A — 1 /r m). Its 10 arcsec long 
slit and 0.30 per pixel sampling in the spatial direction, are 
ideal for observing high-redshift extended astrophysical ob¬ 
jects. The spectrograph was designed to have high through¬ 
put in the blue, with the peak throughput reaching 22 % at 
5600 A including the telescope. 

The quasar was observed in December 2012 with an 
one arcsec width slit aligned along two different position an¬ 
gles (PAs) for 1 hour each. One position was South-North 
(PA = 0) and the other position was East-West (PA = 90). 
Another 1-hour exposure with PA = 90 was taken in Febru¬ 
ary 2013 but the resultant spectrum has a lower SNR. 
Following each exposure, the spectrum of a standard star 
was also recorded allowing us to precisely flux-calibrate the 
quasar spectra. The seeing, measured on the extracted trace 
of the quasar is 1.06 arcsec for PA = 0 and 0.93 arcsec for 
PA = 90. 

We reduced the spectrum using the Mage_Reduce 
pipeline written in the Interactive Data Language (IDL) 
by George BeckeiQ. In addition to the 1-dimensional (ID) 
spectrum, the pipeline provides 2-dimensional (2D) sky- 
subtracted science frames as well as the corresponding 2D 
wavelength solution corrected for the vacuum heliocentric 
velocity shifts. These 2D images will later be used to rectify 
the curved spectral orders (see below). Here, we note that 
since there is an extended Lyman-a emission in order 9 of 
our 2D spectra, we preferentially used a wider extraction 
window to extract this order. 

Two-dimensional spectra when imaged onto a detector 
are often curved with respect to the natural (x,y)-coordinate 
system of the detector defined by the CCD columns (Kelson 
2003). To rectify the curvature of the orders we first rebin 
the wavelength area (given by the 2D wavelength solution) 
and position area (given by the 2D slitgrid array provided 
by the pipeline) increasing the number of pixel by a factor 

^ ftp://ftp.ast.cam.ac.uk/pub/gdb/mage_reduce/ 
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Figure 1. Flux-calibrated spectrum of the quasar SDSS J0823-I-0529 over the full observed wavelength range. The positions of broad 
emission lines are shown with different colors. The Lyman-a-N v broad emission is absorbed by a strong DLA and a strong narrow 
Lyman-o emission is revealed. The cyan curve shows the corresponding error array. The y-axis is in erg s~^ cm~^ A“^. 
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Figure 2. Absorption profiles detected in the DLA at = 3.1910. This redshift is used as the origin of the velocity scale. Multi 
component Voigt profile fits were performed on the absorption profiles and the results are shown by red lines together with individual 
components plotted in different colors. Single component fits were performed on the O I**, Si II*, and Ar I absorption lines. Regions 
corresponding to the A1 II, C II, and O VI are not fitted due to the contaminations by some unidentified absorption lines. 
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of hundred. We define a new 2D array, with one dimension 
representing the wavelength (A) and the other the position 
on the slit (s). For each pixel on the 2D spectrum of the 
quasar, we take its corresponding wavelength and its po¬ 
sition on the slit directly in the rebinned 2D areas. There¬ 
fore, the (x,y) coordinates of each pixel (defined by the CCD 
columns) can now be transformed to an (s. A) coordinate de¬ 
fined by the new 2D array introduced above. In this new (s, 
A)-coordinate system, the curvature of the orders and the 
tilting of the spectral lines are all rectified. In this study, we 
use these rectified 2D images when we discuss the spatial 
extension of the Lyman-a emission line in the quasar spec¬ 
trum. Note that the rebinning process has very little effect 
on the Lyman-a emission because it is conveniently placed 
in the middle of the corresponding order where curvature is 
at minimum. 

Finally, the extracted ID spectrum of each individual 
order is corrected for the relative spectral response of the 
instrument and flux-calibrated using the spectrum of a stan¬ 
dard star (HR 1544) observed during the same night. The 
spectrum of the standard star can be found in the ESO 
standard star catalogue webpag(jf|. We emphasize that our 
standard star spectrum was obtained immediately after ob¬ 
serving the quasar spectrum. We also flux calibrated in the 
same way the 2D spectra in the Lyman-a range. These flux- 
calibrated spectra are then combined weighting each pixel by 
the inverse of its variance. The resulting spectrum (after bin¬ 
ning with a 3 X 3 box) has ~ 27 km s“^ per pixel and 3 pixels 
per resolution element, and therefore FWHM ~ 80 km s“^. 

We have fitted the quasar C iv emission line with two 
Gaussian functions (to mimic the two lines of the doublet) to 
estimate the redshift of the quasar. We derive Zem = 3.1875. 
This is '^330 km s“^ smaller than the redshift of the DLA 
{zdla = 3.1910 derived from the fit of Si ll and Fe ll absorp¬ 
tions) . However, it is well known that redshifts from C iv are 
smaller by up to 600 km s“^ compared to redshifts from nar¬ 
row forbidden lines (e.g. Hewett & Wild 2010). Therefore we 
cannot be certain that the DLA is redshifted compared to 
the quasar. We would need to detect [O ill] lines redshifted 
to the infra-red to have a better estimate of Zem- 


3 A DLA ACTING AS A CORONAGRAPH 

Figure [T] shows the co-added spectrum of SDSS J0823-I-0529 
using all the exposures taken at different PAs. It is apparent 
that a strong DLA is located close to the redshift of the 
quasar. However, strong Lyman-a emission is present in the 
center of the DLA trough. This is possible if the absorbing 
cloud has a transverse extension which is smaller than the 
central narrow line region (NLR) of the quasar. Note that 
the strength of the Lyman-a emission together with the size 
of the cloud derived below makes it very improbable that the 
Lyman-a emission could be a consequence of star formation 
in the DLA itself. 

The redshift of the DLA (^dla = 3.1910) was deter¬ 
mined by conducting a simultaneous Voigt profile fit of the 
Fe II, Si iT, and O T* absorption profiles. 


https://www.eso.org/sci/observing/tools/standards.html 


(b = 10 km s (b = 20 km s 



Figure 3. The relative optical depths of Fe 11 AI 6 O 8 and 1611 
are used to constrain the Doppler parameter of the main DLA 
component. The fit is constrained by Si ll profiles. In each column, 
the Doppler parameter of the strong central component is fixed 
to the value indicated at the top. 


3.1 Elemental abundances 

A Voigt profile fit was conducted on the damped Lyman-a 
absorption line of this system, resulting in a neutral hydro¬ 
gen column density of log Ahi = 21.70 ± 0.10. Note that the 
placement of the quasar continuum is very uncertain as the 
DLA covers the Lyman-a and the N v broad emission lines. 
However, the core and especially the red wing of the DLA 
profile allowed us to satisfactorily constrain the H l column 
density. We will come back to this in Section 3.4. 

We detect absorption lines from O i, O T*, Si ii. Si iT, 
Fe II, Al II, Al III, Ar i, C ii, C iT, C iv. Si iv, P v and 
N V. Absorption profiles are shown in Fig.[2]and results from 
fitting these lines are given in Tables [T]& [2] Techniques used 
here are similar to those in Fathivavsari et al (2013). It must 
be noted that most of the lines are saturated, preventing us 
from deriving accurate column densities especially at the 
resolution of our data (R~4000). 

The profiles are dominated by a main strong component 
clearly seen in particular in Fe 11 AI 6 O 8 and Si ll*A1533. To 
constrain the Doppler parameter of this component, we take 
advantage of the fact that Fe 11 AI 6 O 8 is not saturated and 
well defined while Fe 11 AI 6 II is not detected (see Fig[3]|. We 
start by fitting together Si 11 AI 8 O 8 and 1526, imposing the 
presence of the main component with a fixed Doppler param¬ 
eter. We then use the resulting decomposition to fit the Fe ll 
lines. Results are shown in Fig[3]for 6 = 10 and 20 km/s. It 
is apparent that components narrower than b ^ 20 km/s are 
not favored and b = 10 km/s is definitely rejected. There¬ 
fore, we are confident that imposing 6 ~ 20 km/s for the 
main component in the DLA will give us a good estimate of 
column densities. 

For some of the species, several lines with very differ¬ 
ent oscillator strengths (either in doublets or multiplets) are 
present so that we can derive robust estimates of the column 
densities. This is the case for Si ll, Si iT, and Fe ll. From this, 
we derive metallicities relative to solar, [Si/H] = —0.79 and 
[Fe/H] = —1.87. We also derive an upper limit on N(Fe ill). 
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Table 1. Column densities of low-ionization species. The quoted 
redshift is the redshift of the strongest component in the low ion 
species. 


Redshift 

Ion 

log(N)“ 

Observed 

log(N) 

Model 

3.190974 

Si II 

16.42T0.10 

16.53 


Si II* 

15.49T0.30 

15.49 


Fe II 

15.33T0.20 

16.50 


Fe III 

sS 15.00 

15.09 


Ar I 

14.97T0.20 

15.04 


Al III 

14.80T0.10 

14.64 


O 1“ 

16.45T0.45 

15.70 


O I 

17.08T0.50 

17.62 


“The Doppler parameter is fixed at 20 km s ^ (see Fig. 3). 


We also fit the high-ionization species. We did not try 
to tie the components in different profiles because this high- 
ionization phase could be highly perturbed. In Tabled ex¬ 
cept for the first two components of the C iv absorption 
profiles, all reported column densities are upper limits be¬ 
cause the profiles are strongly saturated. It can be seen how¬ 
ever that the decompositions in sub-components are fairly 
consistent between the different species. 

We detect absorptions from O l** and Si ll* which are 
rarely detected in DLAs (see Noterdaeme et al. 2015, Neele- 
man et al. 2015) and are more commonly seen in DLAs 
associated to GRB afterglows (Vreeswijk et al. 2004; Chen 
et al. 2005; Fynbo et al. 2006). Absorption from the fine 
structure state of Si ll will be used to constrain the density 
of the absorbing cloud (see below). 

The high measured depletion of iron relative to silicon 
(i.e. [Si/Fe] = -1-1.08) in this DLA suggests the presence of 
dust. Consequently, extinction due to dust might be signifi¬ 
cant along this line of sight. Indeed, the median g — r color 
for a sample of 697 non-BAL quasars with redshift within 
/S.Z = ±0.02 around Zem = 3.1910 is {g — r) = 0.30 when 
g — r = 1.1 for QSO J0823±0529. The reddening for this line 
of sight is E(B-V) = 0.09, measured with an SMC reddening 
law template, which places it among the most reddened of 
the sight lines in the Finley et al. (2013) statistical sample. 
We will take this reddening into account in the following 
while discussing the properties of the quasar. 

3.2 Physical conditions in the DLA gas 

In this section we use the photo-ionization code CLOUDY to 
constrain the ionization state of the absorber and its dis¬ 
tance to the quasar central engine. Observed ionic ratios 
of Si II, Si IV, Al III, and Ar l are used to constrain the 
plane parallel models constructed for a range of ionization 
parameters. The calculations were stopped when a neutral 
hydrogen column density of log ^(11 l) = 21.70 is reached. 
Solar relative abundances are assumed and the metallicity 
is taken to be Z = 0.16 Zq from the Silicon abundance of 
the DLA derived in Section 3.1. 

The ionizing spectrum incident on the cloud is taken as 
the combination of the standard ACN spectrum of Math¬ 
ews & Ferland (1987), Haardt-Madau extragalactic spec¬ 
trum (Haardt & Madau 2005, HM05) and the CMB radi¬ 
ation both at 2 = 3.1910. Fig. |4] summarizes the results of 


Table 2. Column densities of high-ionization species. 


Redshift 

Ion 

b 

log(N) 

log(N) 



[km s“l] 

Observed 

Model“ 

3.185037 

C IV 

23 

13.60T0.10 


3.186553 

C IV 

32 

13.90T0.10 


3.189333 

C IV 

86±10 

>15.40±0.10 


3.190953 

C IV 

38±14 

>14.65±0.30 


3.192634 

C IV 

39±6 

>15.40±0.50 


Total(C iv) 



^15.75 

18.43 

3.189333 

Si IV 

56±10 

;sl4.30±0.10 


3.190953 

Si IV 

25±14 

>14.60±0.60 


3.192634 

Si IV 

25±6 

>15.72±0.90 


Total(Si iv) 



^15.75 

16.14 

3.189821 

N v 

99±39 

>14.80±0.30 


3.191157 

N v 

34±29 

^14.89±0.20 


3.192807 

N v 

74±19 

^14.87±0.10 


Total(N v) 



^15.33 

16.55 

3.189788 

P v 

136 

^14.20±0.20 


3.191109 

P V 

97 

>13.95±0.30 


3.192757 

P v 

13 

>13.80±0.60 


Total(P v) 



>14.50 

15.30 


“The Model values are for the ionization parameter logU = —0.3 


Table 3. Hydrogen density (nn) in the cloud, its distance to the 
central AGN (ro), size of the DLA (Ihi) and transverse size of the 
C IV phase for different values of the ionization parameter (U). 


logU 

riH 

[cm-3] 

tq 

[pc] 

m I) 
[pc] 

1{C iv) 
[pc] 

-1.1 

710 

579 

2.3 

3.0 

-0.7 

500 

435 

3.2 

5.8 

-0.3 

355 

326 

4.6 

14.6 

-1-0.0 

250 

275 

6.5 

33.7 

-1-0.3 

180 

229 

9.1 

80.0 


these calculations. Bottom panel gives resulting ionic ratios 
and top panel shows column densities. In both panels mea¬ 
surements and upper limits are indicated by dashed lines. 
For Al II, we scale the Si ll column density assuming solar 
metallicity ratios. 

The column density ratio logAf(Al ll)/A(Al ill) 
yields ionization parameter ranging from —1.1 to 
0.3. Both log A(Si ll)/A(Ar l) ratio and the limit on 
logA(Sill)/A(Silv) are consistent with this range of 
ionization parameters. Our preferred value is log U = —0.3 
and we indicate the column densities for this model in 
Tables 1 and 2. Note that we detect a trough at the position 
of PvA1117 with an absorption profile which is consistent 
with that of other high-ionization species (see Fig. 2). The 
second weaker line of the doublet, PvA1128 is affected by 
noise. The fit to these two lines gives a column density 
which is consistent with the results of the preferred model. 
The same is true for N v. This strongly supports the fact 
that the cloud is highly ionized. 

To determine the hydrogen density of this cloud (i.e. 
jih), a series of CLOUDY models with fixed ionization param- 
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log u 

Figure 4. Results of CLOUDY models of the DLA system. Bot¬ 
tom panel gives resulting ionic ratios and top panel shows colum 
densities. In both panels measurements are shown by dashed hori¬ 
zontal lines and upper limits are indicated with arrows. The thick 
lines in the bottom panel show the range of ionization parameters 
constrained by each ionic ratio. 


eter (varying from log U = —1.1 to 0.3) and varying nn were 
constrncted. Note that in this series of models, the metallic- 
ity and incident radiation are the same as those considered 
above. By increasing the hydrogen density we are indeed 
trying to collisionally populate the excited states of the Si ll 
ground state to explain the observed Si ll* column density. 

Knowing the ionization parameter and the density, we 
can derive the distance of the cloud to the quasar by estimat¬ 
ing the number of ionizing photons emitted by the quasar. 
Since the quasar is reddened, we first estimate the flux at 
20370A (corresponding to H/1 at the redshift of the quasar) 
by extrapolating with a power-law the continuum observed 
at 6125 and 8165 A. Following Srianand & Petitjean (2000), 
we then assume that the de-reddened quasar continuum is a 
power law (/> ~ A“^; with a\ = —1.5) and we consider that 
the flux at ~ 20370 A is not affected by the reddening. We 
thus estimate the flux at the Lyman limit in the rest frame 
of the absorber to be /912A = 3.40x10“^^ erg s“^cm“^A 
This flux corresponds to a luminosity of L912A = 3.25x10'*^ 
erg s“^A“^ at the Lyman limit. We can now estimate the 
rate at which hydrogen ionizing photons are impinging upon 
the face of the cloud by integrating L^/hu over the energy 
range 1 to 20 Ryd. If we assume a flat spectrum for the 



-4000 -2000 0 2000 4000 
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Figure 5. Top panel shows on a velocity scale the quasar C IV 
emission and C iv absorptions from the DLA. The Lyman-a re¬ 
gions for the two PAs (PA = NS in black, PA = EW in red) are 
overplotted in the middle panel. The DLA trough never reaches 
the zero level. Bottom panel shows the subtraction of the two 
spectra (note the change in the y-scale). The residual is consistent 
in the two spectra. Note that the y-axis is in erg s“^ cm“^ A~^. 

quasar over this energy range, we get Q = 6.78x10®® pho¬ 
tons per second. From the definition of the ionization pa¬ 
rameter U, 


dyrr^nHC 

one can see that for given values of nn, U, and Q, one can 
uniquely estimate the distance, ro, from the quasar to the 
absorbing cloud (c is the speed of light). 

We can then derive the size of the cloud along the line of 
sight. For the neutral part, the size along the line of sight is 
I ~ A'hi/u-h. Results are summarized in Table [3] which gives 
the hydrogen density, distance to the quasar and size of the 
neutral phase of the cloud for different values of log U. The 
longitudinal size range from 2.3 to 9.1 pc for a distance of, 
respectively, 579 to 229 pc from the quasar and a hydrogen 
density of 710 to 180 cm“®. Using the structure given by 
the model and assuming that the cloud is spherical we can 
derive the transverse size of the C iv phase. It is given in 
the last column of Tableland range from 3 to 80 pc. These 
estimates for the transverse size of the high-ionization zone 
are only rough estimates because our model is very simple 
(only one density) and we assume spherical geometry. 


3.3 Residual flux in the bottom of the DLA 
trough 

We observe residual flux in the bottom of the DLA trough 
extending in velocity well beyond the narrow emission line. 
This can be seen in Fig. where we overplot zooms of the 
Lyman-a regions observed along the two PAs. It is apparent 
that the flux is never at zero in the trough. We checked in the 
Lyman-a forest that the bottom of saturated lines have on 
an average zero flux. The residual flux is consistent in these 
two spectra as is demonstrated in the bottom panel of the 
figure which gives the difference between the two spectra. 

This means that the DLA cloud does not cover the back¬ 
ground source completely. Since the source of the quasar 
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Figure 6. Upper panel: The residual in the bottom of the DLA 
trough seen over the velocity range [—3000, +3000] km s“^ can be 
explained if ~10% of the broad line region is not covered (solid 
blue line). The de-reddened QSO spectrum (dotted red curve) 
recovered from PCA analysis of the red part of the spectrum is 
used to estimate the quasar emission lines (dashed green line) that 
are then scaled to derive the covering factor of the cloud. Lower 
panel: The red curve is the DLA fit conducted on the de-reddened 
observed spectrum. Green curves indicate the ilcr limits. 


continuum is much smaller than the broad line region (e.g. 
Hainline et al. 2013), this excess is due either to the cloud 
not covering the BLR completely or to a second continuum 
source. 

We investigate here if this excess can be due to an addi¬ 
tional continuum source. If this excess is due to the contin¬ 
uum of the host galaxy then we would expect some residual 
at the bottom of metal lines. Indeed C llA1334 is saturated 
and seems to show some weak residual (see Fig. 2). We thus 
have scaled the LBG mean continuum as given by Kornei et 
al. (2010) to the bottom of C llA1334. We find that this could 
explain 20 to 50% of the DLA residual. The corresponding 
magnitude of the galaxy would be 23.22 which would be 
very bright. In any case this possibility cannot explain all 
the residual. 

It is probable that part of the residual flux in the bot¬ 
tom of the DLA trough is due to the fact that part of the 
Lyman-Q emission from the BLR is not covered by the DLA. 
The size of the broad line region (BLR) can be inferred from 
time delay measurements between variations in the contin¬ 
uum and in the broad lines. Recent investigations of low- 
redshift AGNs show a tight relation between this size and 
the luminosity of the AGN, R = Ax (L/IO^®)®, where R is 
the radius of the BLR, A is a typical distance in light-days 
and L is the luminosity either in an emission line (H/3 or 
G iv) or in the continuum. The index is found to have a 
value close to B~0.6-0.7 when the typical distance A is in 
the range 20-80 light-days for local AGNs (Wu et al. 2004; 
Kaspi et al. 2005). More recently Bentz et al. (2009) find log 
(Rblr) = K + alog(ALA(5100A)) with a = 0.519tQ ggg and 
K = The slope suggests that brighter AGNs have 

to a first approximation the same structure as fainter AGNs 
with only larger dimensions. Therefore, extending the Kaspi 
et al. (2005) relation to higher luminosities yields a typi¬ 
cal radius of the order of 1 pc for the BLR of bright high-z 


quasars. In the present case, luminosity is Lciv ~ 9 x 10^^ 
erg s“^ which gives a size of the BLR of V LI pc. We have 
estimated in the previous Section that the longitudinal size 
of the DLA is of order 2-9 pc. This is larger than the size 
of the BLR. It therefore would mean either that the cloud 
is much smaller in the transverse direction, corresponding 
possibly to a filamentary structure or that some holes are 
present in the cloud or that the cloud is not centered on the 
quasar. 


3.4 QSO BLR Lyman-o: emission and covering 
factor 

To determine the H l column density accurately, one needs 
to reconstruct the shape of the quasar spectrum at the posi¬ 
tion of the DLA profile. We can apply principal component 
(PGA) reconstruction of the quasar flux using the red part 
of the spectrum to estimate the shape of the Lyman-a—N v 
emission (Paris et al. 2011, 2013). To do so, we first subtract 
the residual flux seen at the bottom of the DLA (i.e. resid¬ 
uals from BEL and NEL regions) so that we get zero flux in 
the DLA trough. We then de-redden this spectrum and add 
again the residuals subtracted above. We now have the com¬ 
plete de-reddened observed spectrum. Applying the Princi¬ 
ple Gomponent Analysis (PGA) method on this spectrum, 
we derive the PCA spectrum shown as the dotted red curve 
in the upper panel of Fig. [6] Then, we subtract the contin¬ 
uum from the PCA spectrum (green dashed curve) and scale 
it so that it is consistent with the residual flux seen at the 
bottom of the DLA trough (solid blue curve). This residual 
flux is only ~7 % of that of the Lyman-a broad emission, 
indicating that ~93 % of the Lyman-a BLR is covered by 
the cloud. 

To obtain the flux seen by the DLA, we subtract from 
the dashed green curve the residual flux seen at the bottom 
of the DLA (i.e. the solid blue curve) and re-add the con¬ 
tinuum. We then fit the DLA to obtain the solid red curve 
in the lower panel. We also checked that we get the same 
result if, before fitting the DLA, we first re-redden the PCA 
continuum. As can be seen in Fig. the fit is slightly high 
near the N v emission line although within errors. This is 
possibly because the N v emission line in this quasar may 
be weaker than what is predicted by the PCA method. It 
should be reminded that the PCA reconstruction is an esti¬ 
mate of the quasar spectrum which stays close to the median 
spectrum in the overall quasar population (see Paris et al. 
2011). The result (log A^(H l) = 21.7T0.1) is consistent with 
what was derived previously. 

Figure [7] shows the quasar spectrum as it would be ob¬ 
served if no DLA were present. In this figure, we have added 
the Lyman-a narrow emission component seen in the DLA 
trough to the PCA spectrum. The final spectrum is degraded 
to the SNR ~ 20 by adding Gaussian noise to it. It is appar¬ 
ent from this spectrum that the emission from the narrow 
line region is quite strong in this quasar. 


4 THE QUASAR NARROW LINE REGION 

We call here the narrow line region (NLR), the region of the 
quasar host galaxy that is located within the PSF of the 
observations. This corresponds to about 1 arcsec or 7.1 kpc 
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Figure 7. Reconstructed quasar spectrum as it would be ob¬ 
served if no absorption were present. The original PCA spec¬ 
trum is overplotted as a dotted green curve. The y-axis is in 
erg s“^ cm’ -2A-1. 


at the redshift of the quasar (or a distance of 3.55 kpc on 
both sides of the quasar). The emission seen beyond this will 
be called extended emission. We have shown that the DLA 
is of much smaller dimension so that most of the Lyman-a 
NLR is not covered and is detected as Lyman-a emission in 
the bottom of the DLA trough. We can extract this emission 
which is shown in the top panel of Fig.|8l It can be seen that 
the emission is spread over more than 1200 km s“^ with 
FWHM - 900 km s"A 


4.1 CIV partial coverage 


The large rest equivalent width of the C iv and Si iv ab¬ 
sorption doublets and the flat-bottomed structure of their 
profiles suggest that these lines are saturated (see Fig. [S]). 
However, the flux at the bottom of the C iv doublet ab¬ 
sorption lines apparently does not reach the zero flux level, 
indicating that the absorbing cloud is only partially covering 
the background emission-line region. 

When partial coverage occurs, the residual intensity 
seen at the bottom of absorption lines can be written as 
at each wavelength: 


7(A) = 7o(A)(l - Cf) -b Cf 7o(A) exp [-t(A)] (2) 


where 7o(A) is the incident (unabsorbed) intensity, t(A) is 
the optical depth of the cloud at the considered wavelength, 
and Cf is the fraction of the background emitting region that 
is covered by the absorbing cloud (i.e. the covering factor). 
In the case of doublets, and assuming the covering factor is 
the same for each component of the doublet, we can write: 


1 -|- R2 — 2R2 
1 -b Ri — 2 R 2 


(3) 


where Ri and R 2 are the normalized residual intensities in 
the two absorption lines of the doublet (Petitjean & Sri- 
anand 1999; Srianand & Shankaranarayanan 1999). 

In Fig. [3 we show the covering factor of the C iv. Si iv, 
and N v doublets. Here, the red and blue histograms are 
the profiles of the weaker and stronger transitions of each 
doublet, respectively. The green solid line is the covering 
factor at each point of the profile calculated from Eq.(3). 
The green vertical dotted lines mark the regions used to 
calculate the mean values of Cf, avoiding the wings of the 
profiles. The green horizontal dotted lines along with the 
green filled squares show these mean values. Error bars are 
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Figure 8. Panel a: Lyman-o: emission line profile detected in the 
trough of the DLA. Three Gaussian functions are used to fit the 
profile. Panel b: C iv absorption profiles (the two absorption lines 
C IVA1548 and C ivA1550 are seen). It is apparent from the flat 
bottom of the lines and similar apparent optical depths that the 
covering factor of the C iv gas is not unity. This is demonstrated 
by the model fit overplotted as a red line. Note that a Fe iiA2586 
component at = 1.5140 is blended in the red wing of the 
C IVA1550 profile. Panel c: The residual at the bottom of the 
C IV lines is fitted by scaling the Lyman-a emission profile. The 
red dotted line indicates (on an arbitrary scale) the scaling factor 
(see Eq. (4)). The effective values of the C iv/Lyman-a ratio are 
fitted to be 0.29, 0.045 and 0.18 for the three emission components 
respectively. The C iv/Lyman-a ratio must be much smaller in 
the second emission component. Panel d: Fit of the C iv dou¬ 
blet once the residuals (including both C iv emission lines) are 
removed. Panel e: Same as Panel (c) but the central Lyman-a 
component is supposed to have no C iv associated. Panel f: Fit 
of the C IV doublet once the residuals from panel (e) are removed 
(no C IV is associated with the central Lyman-a component). 


calculated as the standard deviation of Cf in different ve¬ 
locity bins. The measured values of Cf for C iv and Si iv 
are ~ 0.85 and 0.90 respectively suggesting that the size of 
the corresponding BEL region may be similar for the two 
high-ionization species. Note also that N v on the contrary 
seems fully covered. 

The above numbers are however only indicative because 
the resolution of our spectrum is not high (R ^ 4000). In 
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addition there is a possible FeIIA2586 line at Zabs = 1.5140 
blended with the red wing of C lvA1550 (cyan curve in the 
C IVA1550 panel in Fig[2)). In Fig.| 8 ] we demonstrate the ef¬ 
fect of partial coverage on the observed C iv doublet profiles 
taking into account the resolution of the spectrum and the 
presence of the Fell line. The contribution of the Fe llA2586 
absorption to this fit is robustly determined by fitting to¬ 
gether Fe llAA2344,2600 in the same system. In panel of 
Fig. [S] the red curve is the fit conducted with YFFIdlj on 
the original data. Here, we can see that there is no way to 
fit most of the profile without invoking partial coverage. 

4.2 Reconstructing C iv narrow emission line 

(NEL) profile 

As mentioned earlier, the flux at the bottom of the C iv dou¬ 
blet absorption lines does not reach the zero level. This resid¬ 
ual flux could be due to the partial coverage of either the 
C IV BLR or NLR. However, the radius of the C iv phase as¬ 
sociated with the DLA is larger and sometimes much larger 
than the size of the BLR (see Section 3.2) which strongly 
suggests that the residual flux is due to the C IV NLR. 

We thus tried to reproduce the residual by C iv emis¬ 
sion, including both lines of the doublet. Doing this we 
would like to ask the question whether we can associate 
the emission with the Lyman-a emission. Indeed, this could 
bring important clues on the origin of the Lyman-a emis¬ 
sion. In case the Lyman-a emission cannot be associated 
with C IV emission then this would strongly suggest that 
the Lyman-a emission corresponds to scattered light after 
radiative transfer. 

We thus first try to simply scale the Lyman-a profile. 
This is bound to fail as it is apparent that the residuals do 
not follow the shape of the Lyman-a emission. We therefore 
decomposed the Lyman-a profile in three Gaussian functions 
(see top panel of Fig. as suggested by the profile itself. 

This allows us to scale different parts of the Lyman- 
a profile differently. For each pixel we define the ratio 
R — C iv/Lyman-a by combining the three Gaussian func¬ 
tions using the following equation: 

Rl X Crl + i ?2 X Gr 2 -|- Rs X Cr 3 . , 

^ - -Gr + G 2 + G3- 

where Ri, R 2 , and R 3 are the weights for each component 
directly related to the C iv/Lyman-a ratio in this com¬ 
ponent. By assigning different values to Ri, R 2 , and R 3 , 
we can scale each Gaussian individually. For instance, if 
Rl — R 2 — R 3 = 1.0, this results in the red profile shown 
in panel (a) of Fig. [S] which is simply the combination of 
the three Gaussian functions. One can now use the factor 
R = G iv/Lyman-a defined for each wavelength to prop¬ 
erly scale the Lyman-a NEL profile. We thus have to adjust 
the Rl, R 2 , and R 3 parameters until the C iv emission is 
consistent with the residual seen at the bottom of the G iv 
doublet absorption lines. 

We find that the residual can be reproduced (see panel 
(c) of Fig. [S]) with Rl = 0.29, R 2 = 0.045, and R 3 = 0.18. 
The variation of R through the profile is given as a dashed 
pink line in panels (c) and (e) of Fig. [5] The corresponding 

http://www.ast.cam.ac.uk/-rfc/vpfit.html 


fit of C IV after removing the residuals is given in panels (d) 
and (f). 

Since the weight of the second emission component is 
much smaller than the two other ones, we ask the question 
whether it would be possible that the second component has 
no C IV associated. To test this, we impose R 2 = 0.001. The 
result of the fit is given in panels (e) and (f). It is apparent 
that we can find a solution with no C iv associated with the 
second component. 

We thus conclude from all this that (i) the G iv emission 
is strongest around v ~ —200 km s~^; this position could 
indicate the redshift of the quasar, implying that most of the 
Lyman-a emission and the DLA are redshifted; (ii) at least 
part of the NLR Lyman-a emission has no G iv emission 
associated (predominantly around v = -|-100 km s“^) which 
means that the Lyman-a emission in this component is due 
to scattered light or that the emitted gas is located within 
a distance from the quasar smaller than the transverse size 
of the G IV phase associated with the cloud so that the cor¬ 
responding G IV emission is absorbed by the high-ionization 
phase of the DLA cloud. 


5 EXTENDED LYMAN-a EMISSION 

In the middle and bottom panels of Fig. [TU] we show the 
2D spectra of the Lyman-a emission detected in the DLA 
trough for the two PAs. It is apparent that the Lyman- 
a emission is extended and slightly displaced relative to 
the quasar trace. To quantify this we integrate the whole 
Lyman-a profile in the spectral direction and compare the 
result to the spatial PSF derived from the integration of the 
quasar spectrum over the rest of the order beyond 6520 A. 
This is shown in the right panels in Fig. HOI It is apparent 
that the emission is extended well beyond the PSF for both 
PAs and mostly in one direction implying that the total 
emission is shifted towards this direction. 

We investigate whether the extension of the emission 
varies with the velocity position. For this, we split the veloc¬ 
ity range over which the emission is seen in several regions, 
following the profile, and integrate the spatial emission pro¬ 
file over these regions. We then fit the profile by a Gaussian 
function. The results show that along PA = NS the spatial 
extent of the emission is larger than 5 arcsec over about 
2000 km s~^ and that the shift is about 0.2 arcsec towards 
the North direction in the same region. For PA = EW, ex¬ 
tension is about the same but the shift is consistent with 
zero meaning that the extended emission is more symmet¬ 
ric around the trace. The Lyman-a emission is detected up 
to more than 25 kpc from the quasar and there is a strong 
excess emission along PA = NS to the North. 

To better visualize the extended emission, we will sub¬ 
tract the emission associated with the central PSF e.g. the 
emission located on the quasar trace. To do so, the 2D spec¬ 
tral order outside the Lyman-a region of each PA is split 
into several chunks, and in each chunk the counts are in¬ 
tegrated along the spectral axis. A Gaussian is then fitted 
on the profiles to get the central pixel values of the spa¬ 
tial profiles. Finally, we fit a straight line on these values to 
determine the center of the trace for each spectral pixel. 

Once we know the position of the trace exactly, we ex¬ 
tract the spatial profile in each wavelength pixel and we fit 
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Figure 9. Covering factor estimations for C IV, Si IV, and N V doublets following Eq.(3). Vertical dashed lines indicate the windows in 
which it is calculated and horizontal dashed lines (and filled squares) give the mean over these windows. The red and blue histograms 
are the profiles of the weaker and stronger transitions of each doublet, respectively. 


a Gaussian with width equal to the continuum PSF width. 
We then subtract this Gaussian from the profile. Results are 
given in Fig llll Top panels show the resulting 2D spectrum 
of the extension. Bottom panels show the Lyman-a profiles 
integrated along the spectral direction: in red the emission 
which was fitted on top of the trace and in black is the ex¬ 
tension. 

Along PA = NS, most of the extended emission is seen 
to the North up to projected distances of 3.5 arcsec or 
~25 kpc from the center. The striking observation here is 
that the emission is extended over more than 1000 km s-^ 
at all distances from the quasar. Along PA = EW, extension 
is more apparent in one direction as well, to the West. How¬ 
ever the velocity extension gets smaller when the distance 
to the quasar increases. Fig |12l shows the variation of the in¬ 
tegrated flux with distance. In this Figure, the red and blue 
curves indicate variations of flux with distance following the 
relation F = Fo{r/ro)~^ relation. The data roughly follow 
such a law. In case this emission is due to recombination 
and since the ionization flux is decreasing with distance as 
r~^ this would mean that the number of ionized hydrogen 
atoms intersected by the slit is the same at all distances. 
This could be the case if the gas is clumpy so that changes 
in density has little effect. 

The velocity profiles of the extended Lyman-a emission 
at different distances from the AGN are shown in Fig [13] 
and Table [4] summarizes the integrated flux and luminosity 
of the trace and the extension. 


6 DISCUSSION AND CONCLUSIONS 

We have performed slit spectroscopic observations of QSO 
J0823-I-0529 with the MagE spectrograph mounted on the 
Magellan telescope along two PAs, in the North-South and 
East-West directions. The quasar is unique because a DLA 
is located at a redshift very similar to that of the quasar 
{zdla = 3.1910 and 2civ = 3.1875, i.e. Av ~330 km s-^) 
and acts as a coronagraph blocking most of the flux from 
the central regions of the AGN. In the present case, the 
DLA cloud is small enough so that it covers only approxi¬ 
mately 90% of the Lyman-a BLR. This puts us in a unique 
position to be able to directly observe the quasar narrow 
line region and the extended emission line region. Indeed, 


Lyman-a emission is detected up to more than 25 kpc from 
the quasar along both PAs. 


6.1 The quasar NLR 

Finley et al. (2013) have gathered a sample of 57 DLAs 
acting as a coronagraph in front of quasars; the DLAs have 
redshifts within 1500 km s-^ from the quasar redshift. Their 
statistical sample of 31 quasars shows an excess of such 
DLAs compared to what is expected from the distribu¬ 
tion of intervening DLAs. This can be explained if most 
of these DLAs are part of galaxies clustering around the 
quasar. However 25% of such DLAs reveal Lyman-a emis¬ 
sion probably from the quasar host galaxy implying that 
these DLAs have sizes smaller than the quasar emission re¬ 
gion. The Lyman-a luminosities are consistent with those of 
Lyman-a emitters in 75% of the cases and 25% have much 
higher luminosities. The later systems probably reveal the 
central NLR of the quasar. QSO J0823-I-0529 is part of this 
later class of coronagraph DLAs. 

We have shown that the size of the neutral phase is of 
order 2-9 pc which means that the gas does not cover the 
Lyman-a NLR. Three main components of Lyman-a emis¬ 
sion from the NLR are clearly detected at AV = —300, -1-100 
and -1-400 km 3 -^ relative to the DLA redshift (see Fig. 9). 
The corresponding C iv emission seems to be absent in the 
second strongest component. Either the Lyman-a emission 
is scattered or the associated G iv emission is absorbed by 
the high-ionization phase of the DLA. To be covered, the cor¬ 
responding gas must be located within ~3-80 pc from the 
quasar which, we have seen, corresponds to the extension 
of the high-ionization phase of the cloud assuming spherical 
geometry. At least part of the two other components could 
therefore be located beyond this distance from the quasar, 
although part of the third component could also be hidden, 
as the G iv/H l ratio is larger by a factor of about two in the 
first component compared to the third component (see Eig. 
9). Note that the relative velocities of components 1 and 3 
can be interpreted as the result of a conical outflow with 
mean projected velocities of about ±300 km s“^. 
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Figure 10. Panel (a): ID spectrum showing the Lyman-a emission in the bottom of the DLA trough. Black and red histograms are for 
PA = NS and EW, respectively. Red, blue and green horizontal dotted lines mark the zero flux level, 1 cr, and 3 cr uncertainty on the 
flux, respectively. Panels (b) and (c): 2D spectra of the Lyman-a emission seen in the trough of the DLA at ^abs = 3.1910 for PA = NS 
(panel b) and PA = EW (panel c). The contours show the observed flux density with a separation between two lines of 3.5 x 10“^^ 
erg cm~^ A The outermost contour corresponds to a flux density of 4.37 x 10“^® erg s“^ cm“^ A The black horizontal 
dotted line shows the center of the continuum trace. The blue and red dashed horizontal lines show the FWHM of the trace and extended 
Lyman-a emission, respectively. Panels on the right show the spatial emission profile integrated over the whole wavelength range (blue 
curves) together with the continuum PSF (red lines) derived from the red part of the order. Note that in Panel (a), the y-axis is in 
erg s“^ cm”^ A 
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Figure 11. Upper panels: 2D spectra of the Lyman-a emission seen in the trough of the DLA (left PA = NS and right PA = EW) 
after subtracting the emission in the PSF centered on the quasar trace. The contour levels are the same as in Fig llOl The red horizontal 
dotted line shows the center of the trace and the red solid horizontal lines show the FWHM of the trace. Bottom panels: Emission profiles 
integrated over the whole wavelength window of the Lyman-a emission in the PSF centered on the quasar trace (red histograms) and the 
extended emission (black histograms). Note that in the upper (resp. lower) panels, the y-axis is in arcsec (resp. erg s“^ cm~^ arcsec”^). 


6.2 Lyman-a emission from the center to the 
outskirts of the host galaxy 

Extended Lyman-a emission has been observed around 
high-redshift radio galaxies (Heckman et al. 1991; van Ojik 
et al. 1996), as well as around radio-quiet quasars (Bunker et 


al. 2003; Christensen et al. 2006). For radio-loud quasars or 
radio-galaxies, the Lyman-a flux of the nebula is an order of 
magnitude higher (Christensen et al. 2006), presumably be¬ 
cause the emission of radio-loud quasar gaseous envelopes is 
enhanced by interactions with the radio jets. More recently, 
North et al. (2012) have used a careful treatment of the spec- 































































12 H. Fathivavsari et al. 


Table 4. Lyman-o integrated fluxes and luminosities in the PSF 
centered on the quasar trace and in the extension. The last row 
gives the integrated flux and luminosity of the broad Lyman-a 
emission calculated using the simulated spectrum shown in Fig[7] 



Integrated flux 

Luminosity 


[erg cm“^] 

[erg s-l] 

Trace(PA=NS) 

3.59 X 10-1® 

3.43 X lO"!® 

Extension (PA=NS) 

3.97 X 10“i® 

3.79 X lO'i® 

Trace(PA=EW) 

3.36 X 10-1® 

3.21 X 10"!® 

Extension(PA=EW) 

2.68 X 10-1® 

2.56 X 10"!® 

Lyman-a (BELR) 

7.10 X 10-1® 

6.78 X lOii 


tral PSF to extract quasar traces. This revealed four detec¬ 
tions of extended emission out of six radio-quiet quasars at 
z ~ 4.5 with extensions of diameter 16 < d < 64 kpc down 
to a luminosity of 2x 10“^^ erg s“^ cm“^ arcsec”^. The emis¬ 
sion has 900 < FWHM < 2200 km s“^. Our observations 
are in line with these numbers. The extended emission we 
detect in QSO J0823-I-0529 has a diameter of ~50 kpc and 
FWHM ~ 900 km s“^. However, QSO J0823+0529 does 
not seem to follow the L(Lyman-a) vs L(BLR) relation in¬ 
dicated by these authors. Indeed the Lyman-a luminosity in 
QSO J0823-I-0529 is more than an order of magnitude larger 
than what would be expected from this relation even after 
correcting for dust attenuation. It is still possible that dust 
Is present closer to the quasar and further attenuates the 
BLR Lyman-a emission. It is also possible that the decom¬ 
position between the broad and narrow line region emissions 
was ambiguous in previous studies so that the NLR emission 
could have been underestimated. 

The emission is more extended to the North-West of 
the object (see Figs. 11, 13 and 15). There are two notable 
features in the spatial and velocity structure of the neb¬ 
ula. First, in the North-West direction, the kinematics are 
strikingly similar along the trace and 10 kpc away from the 
center (see Fig. 13) with a velocity spread of more than 
1000 km s“^. The emission is quite strong in this region. 
Secondly, there Is a clear gradient of about 1000 km s~^ 
between 15 kpc to the East and 20 kpc to the West. It is 
tempting to interpret these features as the superposition of 
the emission of gas in the disk of the galaxy, where the den¬ 
sity Is higher and turbulent kinematics prevent gas clouds to 
be well organized with emission from gas flowing out of the 
disk with velocities of the order of 500 km s“^. This gas is 
best seen up to 20 kpc to the West and 10 kpc to the East. 
Such winds can be reproduced by recent models taking into 
account the effects of radiation trapping (Ishibashi & Fabian 
2015). 

6.3 Nature of the DLA 

It is well demonstrated that bright quasars are surrounded 
by large amounts of gas both in extended ionized halos up 
to 10 Mpc from the quasar (e.g. Rollinde et al. 2005) and 
extended (300 kpc) reservoirs of neutral gas in the halo of 
the host galaxy (Prochaska et al. 2013, 2014; Johnson et 
al. 2015). However, along the line of sight to quasars, the 
incidence of neutral gas is less (e.g. Shen & Menard 2012) 
indicating that the ionizing emission from quasars is highly 
anisotropic. The DLA we discuss here is part of the so-called 



Position (kpc) 



Figure 12. Variations of the integrated Lyman-o flux with dis¬ 
tance from the AGN. The black lines with error bars are the 
observations and the red and blue curves indicate a variation of 
the flux with distance following a r~^ relation. Note that in both 
panels Fq is in erg s“^ cm“^ arcsec“^. 


proximate DLAs in the sense that the absorption redshift is 
similar to that of the quasar (zabs ~ Zem, see e.g. Ellison et al. 
2010). However, this is the first time it can be demonstrated 
that the gas is located very close to the quasar (< 400 pc) 
and probably associated with the central part of the host 
galaxy. 

We have modelled the physical state of the gas in the 
DLA using photo-ionization models. The ionization param¬ 
eter is found to be in the range —1.1 < log!/ < -1-0.3 which 
means that the cloud is mostly highly ionized. A density of 
riH ~ 710-180 cm“® is needed to explain the absorption lines 
from O I and Si ll ground state excited levels implying that 
the neutral phase should have dimensions of the order of 2 
to 9 pc and be embedded in an ionized cloud of size ~3- 
80 pc. From this, we could derive that the cloud Is located 
between 230 and 580 pc from the quasar. The metallicity 
of the cloud, Z = 0.16 Zq, is typical of the metallicity of 
standard intervening DLAs (Rafelski et al. 2012). 

It is intriguing to note that the high-ionization phase 
of the cloud we see has characteristics similar to those of 
some of the warm absorbers seen in many AGNs. Tombesi 
et al. (2013) argued that warm absorbers (WA) and ultra¬ 
fast outflows (UFO) could represent parts of a single large- 
scale stratified outflow observed at different locations from 
the black hole. The UFOs are likely launched from the in¬ 
ner accretion disc and the WAs at larger distances, such as 
the outer disc and/or torus. There are still significant uncer¬ 
tainties on the exact location of this material, which ranges 
from a few pc up to kpc scales, (e.g. Krolik & Krlss 2001; 
Blustln et al. 2005). The absorption lines are systematically 
blue-shifted, indicating outflow velocities of the WAs in the 
range 100-1000 km s“^. King & Pounds (2013) have argued 
that the dense gas which surrounds the AGN when it starts 
shining is swept out by the fast winds powered by the ac- 
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cretion luminosity. The wind is halted by collisions near the 
radius where radiation pressure drops. The shocked gas must 
rapidly cool and mix with the swept-up ISM. Distance from 
the AGN and properties of the gas are similar to what is 
observed for warm-absorbers, as in our case. Therefore, in 
QSO J0843-I-0529, the DLA could be located in the galac¬ 
tic disk at the terminal position of the wind, at the limit of 
the interstellar medium of the host galaxy. The presence of 
dust in the gas can be considered as supporting this view as 
it is expected in the dense environment of AGNs (see e.g. 
Leighly et al. 2015). 

The only caveat with this idea is that the metallicity 
of the DLA is typical of intervening DLAs, when we would 
expect the gas in the ISM of the quasar host galaxy to have 
larger metallicity. This, together with the low outflow veloc¬ 
ity (the DLA is centered on the Lyman-a emission), argues 
for another explanation for the origin of this gas. Given the 
distance to the central AGN (230 < tq < 580 pc), one is 
tempted to conjecture that it could be the final fate of in¬ 
falling gas. Note that in that case, the presence of dust is 
not a problem as the observed amount is not unusual com¬ 
pared to what is seen in typical high column density DLAs. 
Indeed, accretion is believed to happen through cold flows 
which are expected to be of lower metallicity compared to 
the environment of the quasar (see Bouche et al. 2013). Lit¬ 
tle is known about how quasars at high redshift interact 
with cold infalling streams of gas, and in particular whether 
these collimated structures can survive the energy released 
by the AGN. Therefore the DLA cloud could be in a tran¬ 
sitory phase before it is completely destroyed by the AGN 
ionizing radiation field (Namekata et al. 2014). 

The two alternatives are actually not incompatible, 
since some of the warm absorbers could be part of the galac¬ 
tic ISM swept away from the center by the AGN and placed 
at large distances from the AGN (e.g. King & Pounds 2013). 
These outflowing warm absorbers may then intersect with 
accreted cold streams. 


6.4 Concluding remark 

There is little doubt that DLAs acting as coronagraphs are 
important targets to be observed and analyzed. Besides re¬ 
vealing interesting characteristics of the quasar host galaxy 
they are, like in QSO J0823+0529, potentially part of the 
machinery that power the AGN. Finley et al. (2013) found 
57 such systems out of about a third of the BOSS targets. 
With the on-going eBOSS-SDSSIV survey and the forseen 
DESI survey (Schlegel et al. 2011), the number of such sys¬ 
tems will increase by a large factor in the future. It will thus 
be important to gather a large sample of these objects to be 
able to study their characteristics statistically. 
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Figure 13. Velocity profiles of the extended Lyman-a emission at 
PA = EW (right panels). The parameter rm gives the mean distance 
shows the spectrum of the trace for the PA = EW (resp. PA = NS). 



-2000 -1000 0 1000 2000 
Velocity (km s~^) 


different distances from the AGN for PA = NS (left panels) and 
to the AGN. The red histogram in the left panel (resp. right panel) 
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